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ABSTRACT 

The Atmospheric Imaging Assembly onboard the Solar Dynamics Observa- 
tory observed a coronal mass ejection with an embedded filament on 2011 Febru- 
ary 24, reavealing quasi-periodic vortex-like structures at the northern side of 
the filament boundary with a wavelength of approximately 14.4 Mm and a prop- 
agation speed of about 310 ± 20 km s _1 . These structures could result from 
the Kelvin-Helmholtz instability occurring on the boundary. We perform 2.5D 
numerical simulations of the Kelvin-Helmholtz instability and compare the sim- 
ulated characteristic properties of the instability with the observations, where 
we obtain qualitative as well as quantitative accordance. We study the absence 
of Kelvin-Helmholtz vortex-like structures on the southern side of the filament 
boundary and find that a magnetic field component parallel to the boundary 
with a strength of about 20% of the total magnetic field has stabilizing effects 
resulting in an asymmetric development of the instability. 

Subject headings: Instabilities — Magnetohydrodynamics (MHD) — Methods: 
numerical — Sun: corona — Sun: coronal mass ejections (CMEs) 



Introduction 



For the dynamics of the solar atmos phere, magneto 



are thought to play a ma j or role 



Ofman fc Thompson 



2011 



Taroyan 



e^ 



Ryutova et a' 



2011 



Innes et al. 



lydrodynamic (MHD) ins tabilities 



20101 : 



Foullon et al. 



2011 



20121 ). Specifically, there are 



boundary layers, where the Kelvin-Helmholtz instabilit y (KHI) might be able to deve 



treated in several theoretical studies on this top i c (e.g. 



Zaqarashvili et al. 



20101 : iSoler et al 



Karpen et al.lll993l: lAndries and Goossens 
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2001 



Hevvaerts fc Priest 



Kolesnikov et al 



2004 : 



1983 



Rae 



Soler et al. 



op, as 



1983 



20101 : 



20121 ). In principal, the KHI a rises at boundaries with 



96 111, and t here exist many 



Miura 



19971 . and references 



a velocity shear parallel to a boundary layer ( jChandrasekhar 
plasma configurations in space where the KHI can occur (e.g. 
therein) . 

With the unp recedented observations by the Atmospheric Im aging Assembly (A IA; 



Lemen et al. 



Pesnell et al. 



20121 ) 



20121 ) onboard Solar Dynamics Observatory (SDO; 
it is possible to gain observational evidence for instabilities in the solar atmosphere. 
SDO/AIA carries four telescopes providing narrow-band full-disk imaging of the Sun in 



different UV and 



1".5 ( ILemen et al. 



UV p assbands at a cadence as h igh as 12 s and spat i al reso lution of 



20121). 



Foullon et al 



( I2OIII ) and lOfman fc Thompson! ( J201 lh were the 



first to present obser v ation s of the KHI associated with a coronal mass ejection (CME). 



Ofman fc Thompson! ( 1201 lh analyzed a series of vortices seen along the boundary of an 
evolving coronal dimming region on 2010 April 8. They compared for the first time the 
observational features with the result s of a 2.5D MHD model of the KHI and found good 



qualitative agreement. 



Foullon et al. 



( j201ll ) reported observations of wavy vortex-like 
structures on one flank of a CME that occurred on 2010 November 3. These studies hint 
at the existence of the KHI in the solar corona. Yet, we need to have more observations 



and theoretical (numerical) studies on this subject to be clearly able to distinguish between 



structures produced by the KHI and such that have a different origin. 



In this letter, we report quasi-periodic vortex-like structures observed on the boundary 
of an erupting filament /CME on 2011 February 24 in SD O/AIA high cadence imagery. We 
follow a similar approach as lOfman fc Thompson! ( 1201 ll ) and give a detailed comparison 
of these observations and a 2.5D MHD simulation of the KHI. Qualitative as well as 
quantitative analyses of both spatial and temporal features are presented. Furthermore, 
we investigate the effect of a magnetic field component parallel to the boundary layer on 
the evolution of the instability and determine a threshold for the field strength needed to 
stabilize the boundary. 



2. Observations and Analyses 



On 2011 February 24, a CME eruption with an embedded filament was observed 
by the AIA instrument onboard SDO in association with a GOES M3.5 limb flare 



( jBattaglia and Kontar 



2011 



Martinez Oliveros et al. 



20121 ) . The top three panels of Figure 



CD show direct images in the AIA 304 A passband (dominated by emission from He II; 
log (T) = 4.7), the bottom panels display running difference images i n the 171 A channel 



Lemen et al. 



20121 ). On the right panel 



(dominated by emission from Fe IX; log (T) = 5.8; 
we show the associated CME observed by the Large Angle and Spectrometric Coronagraph 
(LASCO) C2 onboard the Solar and Heliospheric Observatory (SOHO). The insert (yellow 
rectangle) represents the AIA 171 A running difference image recorded about 30 minutes 
before the LASCO image shown. The CME structure breaks up rather early in its evolution 
(around 07:30:00 UT), and quasi-periodic vortex-like structures develop around 07:36:00 
UT on the northern side of the boundary of the filament embedded in the CME (see 



accompanying movie in the online version of this publication) . 
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Fig. 1. — Left: Observations of SDO/AIA of the erupting filament. The 304 A channel 
shows direct images (movie in online version), the 171 A channel the corresponding running 
difference images. Right: Observations of SOHO/LASCO of the associated CME. The 
yellow rectangle shows the AIA 171 A running difference image about 30 minutes prior to 
the LASCO observations. 
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Fig. 2. — Zoom into the AIA 304 A images, showing the northern boundary of the filament 
and indicating the height of the observed structures. 



The vortex structures are best observed in the AIA 304 A passband (T ^ 0.5 MK, see 
Figures [1] and Ej) . Wavebands imaging the outer corona in temperatures between 0.4 MK 
- 2 MK (131 A, 171 A, 193 A and 211 A) show only the parts of the vortices close to the 
boundary layer. 

Measuring the front of the filament, we derive an average speed of w 470 km s _1 . The 
height of the structures lies in the range of 2.5 to 4 Mm (Figure EJ). In Figure [3] we show 
cuts along the northern edge of the filament for different times. These plots clearly reveal 
the quasi-periodic appearence of the vortex-like structures, their evolution and propagation 
with time, which allows us to estimate characteristic quantities from these observations. 
The spatial period of the vortex-like structures is found to be approximately 20" (Figure EJ), 
which corresponds to a wavelength A of approximately 14.4 Mm. The average propagation 
speed ^ prop of the observed structures along the boundary is about 310 ± 20 km s _1 . Since 
this event occurred on the limb of the Sun, we can neglect projection effects. 

3. Simulation Setup 



We solve the ideal M' 



scheme (Toth fc Odstrcil 



3D equations numerically with the so-called TVD Lax-Friedrichs 
1999 ) . The numeri cal algorithm was already used in previou s 



studies on the Kelvi n-Helmholtz instability ( jAmerstorfer et al. 



Zellinger et al. 



2012 ). 



2010 : 



Mostl et al. 



2011 



We estimate the pl asma parameters in the corona using typical values, which can be 



found in the literature (lAschwanden 



20041 ) . In the corona, the electron density n e is taken 
to be 1 x 10 15 m -3 , the plasma pressure p is 0.09 Pa, and the magnetic field B is 10 x 10 -4 
T. With these values, we obtain for the total pressure II = 0.5 Pa, which is the sum of 
plasma and magnetic pressures. In the filament, we assume that the density is at least 10 
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Fig. 3. — Profiles of the quasi-periodic structures along the northern boundary layer of the 
filament for different times, measured from direct images in the AIA 304^4 channel. The 
time is given in UT. The propagation speed of these structures is about 310 ± 20 km s _1 . 



times higher than the corona value and that T < 1 x 10 5 K. With further assumptions that 
rii « n e = n and T^ T e = T we have for the plasma pressure p = 2nkT = 0.028 Pa. The 
subscript i denotes the values for the ions. We take 450 km s _1 as the approximate speed 
of the filament Vf and assume that the surrounding coronal plasma moves at half of this 
speed. These plasma parameters are used as input for the simulations. 

The general setup for the simulations is such that the x-direction is along the boundary 
layer and the ^/-direction is perpendicular to it. This means that the plasma parameters 
change in y-direction, and the disturbances (vortices) propagate along the x-direction. The 
background plasma parameters vary according to a tanh-profile. Thus, we do not have a 
tangential discontinuity, but the boundary layer has some finite thickness. 

As initial condition we need an equilibrium plasma configuration. In ideal MHD 
simulations, this means that the total pressure II has to be constant across the boundary 
layer, i.e. in our case 

with 

n = P + ^, (2) 

where p is the thermal pressure, B is the magnetic field and {1$ is the permeability in 
vacuum. In order to have an initial equilibrium, we determine the total and the plasma 
pressure from the plasma background parameters and then calculate the magnetic field 
according to 



B = y/2fJL (U- P ). (3) 

For the cases, where we have two magnetic field components, B x and B Z) we assume a 
constant B x and then use Equation [3] to calculate B z . Onto this equilibrium we impose 
small random initial disturbances of v x and v y , with an amplitude of 1% of the initial 
filament velocity. 



We use normalized quantities in the numerical simulations. The normalization 
constants are the initial filament speed Vf ) the filament mass density pf and the half- width 
of the boundary layer a. Accordingly, the time is then normalized with a/vf, the pressures 



with Pfvj and the magnetic field with J ji^pfV 2 ,. The computational box extends from 
to 50 a in x-direction and from —25 a to 25 a in ^-direction. At the x-boundaries we impose 
periodic boundary conditions and at the ^/-boundaries transmissive ones. 

4. Results 

The top panel of Figured] shows the (normalized) mass density for different (normalized) 
times in the simulation run. The blue layer represents the surrounding coronal plasma, 
whereas the red layer is the ejected filament. For this simulation run, only a z-component of 
the magnetic field was taken on both boundaries. The x- and y-components of B are zero, 
as well as the z-component of the velocity. The plasma flows from right to left. After some 
time, small initial perturbations grow in amplitude and eventually reach their nonlinear 
vortex phase. Then, a saturation sets in and the regular vortex-structures become irregular, 
as does the whole boundary layer. At t w 200 the vortices have fully formed on both 
boundary layers. However, the observations show vortices only on one side of the filament 
boundary. 

From theory it is known that a magnetic field parallel to the boundary layer can have 
a stabilizing effect on the KHI due to magnetic field line tension. The Alfven Mach number 
Ma ) determined only with the parallel magnetic field component, has to be larger than 2 
for the instability to develop. Combining the effects of compressibility and magnetic field 
line tension results in the instability condition that 2va < Av < 2v s , where Av is the 
total velocity jump across the boundary layer, va is the Alfven velocity and v s the sound 
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Fig. 4. — Time series of the mass density (color code). The plasma flows from right to 
left. Top panel: Only a z-component of the magnetic field was assumed. Middle panel: An 
x-component of the magnetic field of 10% of the total magnetic field strength was assumed. 
Bottom panel: An x-component of the magnetic field of 20% of the total magnetic field 
strength was assumed (movies in online version). 



ve. 
to 



ocity, or equivalently 



Miura and Pritchett 



11 



va < v s (for the exact derivation of the instability criterion we refer 



19821 ) . We performed simulations including an x-component of the 
magnetic field at the lower boundary layer to derive the field strength needed to stabilize 
the boundary layer. 

The middle and bottom panels of Figure @] show the results of two simulation runs 
where an x-component of the magnetic field is included at the lower boundary layer. The 
middle panel is the result with B x = 0.1 B and the bottom panel with B x = 0.2 B. For the 
first case, there are still vortices visible at the lower boundary, whereas for the latter case, 
the parallel magnetic field is large enough to stab ilize the boundary layer. This stabilizing 



effect of a strong parallel B ( jChandrasekhar 



1961 



is a possible reason why the vortices are 



only observed on one side of the erupting filament. 

From the simulation that most closely resembles the observations (bottom panel of 
Figure Hj), we determine the wavelength A (i.e. the distance between two vortices) and the 
height h of the vortex-structure. We derive an approximate A of about 10 to 15 a and an 
approximate h of about 3 to 5 a. 

In order to get physical values, we have to make an estimation of the half width of the 
shear layer a. From the observations shown in Figure El we can estimate the boundary layer 
to be about 1 to 2 Mm, resulting in a w 0.5 — 1 Mm. Taking these values, the wavelength 
and the height of the simulated structures yield (i) A ~ 5 — 7.5 Mm and h ~ 1.5 — 2.5 Mm, 
or (ii) A w 10 — 15 Mm and h w 3 — 5 Mm, respectively. The propagation velocity ^ pr0 p of 
the simulated KH vortices lies around 330 km s _1 . All of these characteristic features of the 
simulated instability are in basic quantitative agreement with the observations (A w 14.4 
Mm, h w 2.5 — 4 Mm, ^ prop w 310 ± 20 km s _1 ). Furthermore, we find that a B x of 
approximately 2 x 10 -4 T is needed to stabilize the lower boundary layer. 

Qualitatively, we can compare the structure of the simulated vortices with the 



12 



observations. Figure [5] shows on the left side the evolution of the simulated KH vortices and 
on the right side the profiles through the vortices for each time step. A rather characteristic 
double peak appears when the vortices are rolled-up. Such double peaks are also visible 
in the profiles obtained from the direct images presented in Figure [31 Another feature 
associated with an instability is a growing amplitude of the perturbations, which is seen in 
both, simulations and observations (Figures [5] and El respectively). 

5. Discussion and Conclusions 



There are various impor tant implications res ulting from the occurrence of the KHI at 



the boundaries of CMEs. As 



affect the drag force ( jVrsnak et al 



Foul 



on et 



al. 



( 1201 lh discuss, the KHI at CME surfaces may 



2012 ). and hence, the CME kinematics in interplanetary 



space and its geoeffectiveness. Due to rec onnection inside KH vo rtices, a plasma transport 



across the boundary can be initiated (e.g. iNykyri and Otto 



2004 ) leading to a mass loss of 



the CME. Additionally, magnetic reconnection can le ad to the erosion o f the initial magnetic 



flux, again influencing the g eoeffectiveness of CMEs (jDasso et al. 



20ld; 



Ruffenach et al. 



20121) 



2007 : 



Taubenschuss et al. 



We present numerical simulations of the KHI on a boundary of a filament embedded 
in a CME. We find qualitative as well as quantitative agreements between the simulation 
results and the observed features of quasi-periodic vortex-like structures observed at the 
northern boundary of a filament during an eruption on 2011 February 24. Additionally, we 
have shown that a magnetic field component parallel to the boundary layer with a strength 
of about 0.2 S, corresponding in our case to 2 x 10 -4 T, results in a stabilization of the 
boundary layer and could be an explanation of the observed asymmetric development of the 
vortex-like structures. Our assumption of a tenfold density jump towards the filament is a 
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Fig. 5. — Left: Evolution of the KH vortices. The black lines indicate, where the cuts shown 
on the right plots were taken. Right: Profiles of the density through the vortices at different 
time steps. 
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very conservative one, as one c an expect a larger 



was shown in previous studies ( jAmerstorfer et al. 



density increase o: 



2010; 



Mostl et al. 



about 100 or more. As 



20111 ). a larger density 



ratio results in smaller instability growth rates and, thus, it exhibits a stabilization of the 
boundary layer. Contrary to this stabilizing effect, a larger velocity shear than the one 
assumed in this case study would act destabilizing. 

The numerical MHD model applied is of course an idealized one with simplifying 
assumptions of the geometry of the CME and of the parameters. The model uses standard 
input plasma parameters, which are typical for the observed features but may differ for this 
specific event. Nevertheless, it captures important basic physical properties of the studied 
boundary layer processes. We point out that we are limited in observations coming from 
a single vantage point. The vortex structures might be affected by projection effects and 
possibly signatures of twisted field lines (flux rope), viewed from aside, may show similar 
structures. 

The stabilizing effects might be large enough, such that the KHI is only seldomly 
able to develop at CME boundaries in the solar corona. For the present case study we 
found evidence for the occurrence of the KHI on boundary layers in the solar corona. More 
observational data of possible KHI formation are needed to better understand its physics 
and its consequences on propagating MHD structures. 

UVM, MT and AMV gratefully acknowledge the Austrian Science Fund (FWF): 
P21051-N16, V195-N16 and P24092-N16. 

Facilities: SDO (AIA) 
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